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Abstract 

Oh 



We construct an effective non-relativistic quantum field theory that describes 
bound states of vr+vr" pairs and their hadronic decays. We then derive a gen- 
eral expression for the lifetime of the ground state at next-to-leading order 
in isospin breaking. Chiral perturbation theory allows one to relate the de- 
5^ \ cay rate to the two S'-wave vrvr scattering lengths and to several low-energy 

constants that occur in the chiral Lagrangian. Recent predictions for the 
scattering lengths give r = (2.9 ± 0.1) x 10~^^ sec. This result may be con- 
fronted with vr^vr" lifetime measurements, like the one presently carried out 
at CERN. 
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I. INTRODUCTION 



The DIRAC experiment at CERN |^ aims to measure the hfetime of the vr+vr^ atom 
(pionium) in its ground state with high precision. This atom decays predominantly into 
two neutral pions, F ^ ^271^' ■ The latter decay rate is proportional to the square of the 
difference oq — CI2 of the strong S-wave vrvr scattering lengths with isospin J = 0, 2. 
The measurement will therefore allow one to determine this difference, which may then 
be confronted with the predicted value ao — 02 = 0.265 ± 0.004 [Q. What makes this 
enterprise particularly exciting is the fact that one may determine in this manner the nature 
of spontaneous chiral symmetry breaking in QCD by experiment: Should it turn out that 
the predictions are in conflict with the results of DIRAC, one would have to conclude 
that spontaneous chiral symmetry breaking in QCD differs from the standard picture 
An analogous determination of the nature of spontaneous chiral symmetry breaking may be 
performed through an analysis of K^,^ decays, which allows one to measure the scattering 
length ao 0. 

In order to determine the scattering lengths through a measurement of the pionium life- 
time, the theoretical expression for the width must be known with a precision that matches 
the accuracy of the lifetime measurement of DIRAC. In Ref. (see also we have 

presented a compact expression for r27rO in the framework of QCD (including photons) by 
use of effective field theory techniques. The result obtained contains all terms at leading 
and next-to-leading order in the isospin breaking parameters a ^ 1/137 and (m„ — m^)^. 
On the basis of this formula, a numerical analysis was carried out in Ref. [|12| at order e^p^ 
in ChPT. The aim of the present paper is: i) to give a complete description of the theory of 
the TT+TT^ atom decay, providing the details that were omitted in [ll0| , p!2| , and ii) to update 
that numerical analysis by use of the information recently obtained in ||^ on the scattering 
lengths and on one of the low-energy constants. 

We first briefly review previous work on the subject. Theoretical investigations of 
hadronic atoms and, in particular, of 7r~^n~ decay, have been performed in several set- 
tings. Potential scattering theory in the framework of quantum mechanics has been used 
in PJTB|-P^, and methods of quantum field theory have been invoked as well [|T^-|2T[]. In par- 
ticular, in Refs. 0], the lifetime of the tt+tt^ atom was calculated by use of two-body wave 
equations of 3D-constraint field theory. In Refs. |2T[], the 7r"'"7r~ atom decay was studied in a 
field-theoretical approach based on the Bethe-Salpeter equation. The results for the 7r"'"7r~ 
atom lifetime obtained with the two latter approaches contain the major next-to-leading 
order terms in isospin breaking and agree both conceptually and numerically. However, in 
these investigations the momentum dependence of the strong tttt scattering amplitude was 
neglected. 

In several recent publications |[2^-p7|, the decay of tt+tt^ atoms has been studied in 
the framework of a non-relativistic effective Lagrangian - a method originally proposed by 
Caswell and Lepage |2^ to investigate bound states in general. This method has proven to 
be far more efficient than conventional approaches based on relativistic bound-state equa- 
tions. It allows one e.g. to go beyond the approximation used in pO| , pl| for the scattering 
amplitudes. In our previous publications [0,|ll], we have used the same method. We refer 
the reader to for a comparison of the various results obtained in the effective framework. 

We now describe the general features of the system that we are going to study. The 
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vr+vr" atom is a highly non-relativistic, loosely bound system. The pions are mainly bound 
by the Coulomb force, and the atom decays predominantly through the strong interactions. 
The average momentum of the constituents in the CM frame is ~ 0.5 MeV, and the Bohr 
radius of the bound state is ~ 400 fm. The decay width of the tt+tt" atom ~ 0.2 eV is 
much smaller then the binding energy ~ 2 x 10^ eV. For this reason, a non-relativistic 
framework provides the most economical and powerful approach to the calculation of the 
characteristics of this sort of bound states. Since the strong interactions between pions at 
low energy can be described with ChPT, the theory of the 7r"^7r~ atom turns out to be a 
merger of a non-relativistic approach with ChPT. Owing to the might of the non-relativistic 
approach which almost trivializes the calculations in the bound-state sector, we are able to 
determine the first few coefficients in the chiral expansion of the bound-state observables. 

The paper is organized as follows. In Section |T| we discuss the foundations of the theory: 
the non-relativistic Lagrangian, Green functions, and matching to the relativistic ampli- 
tudes. Bound states are discussed in Section |TTI|. Using Feshbach's formalism we derive 
a master equation for the position of the poles in the resolvent. In Section |I^ we derive - 
on the basis of the master equation - a general expression for the decay width of the tt+tt" 
atom in the ground state, valid at next-to-leading order in isospin breaking. We then express 
this quantity, through the matching condition, by the relativistic 7r"'"7r~ — > vr^vr^ scattering 
amplitude at threshold. A numerical analysis of the decay width at order e^p^ in ChPT is 
also carried out in this section. Section ^ contains our conclusions. Background material is 
relegated to the Appendixes: In Appendix 0, we discuss the construction of a general non- 
relativistic Lagrangian with pions and photons. The scattering sector of the non-relativistic 
theory is discussed in Appendix |B[ In particular, we argue that the contributions of trans- 
verse photons to the tt+tt" vr^vr^ scattering amplitude vanish at threshold at order e^. 
Therefore, these diagrams may be omitted in matching the relativistic and non-relativistic 
amplitudes. Appendix P deals with the bound states in the non-relativistic theory: we 
show that, for a large class of diagrams, transverse photons do not contribute to the decay 
width at next-to-leading order in isospin breaking. On the basis of the results obtained 
in Appendixes ^ and P, we completely eliminate transverse photons from the theory. In 
Appendix we compare two different matching procedures. Finally, in Appendix ^ the 
SU (3) X SU (3) SU (2) X SU (2) mapping of the pertinent combination of electromagnetic 
low-energy constants in ChPT is provided. 



II. THE EFFECTIVE NON-RELATIVISTIC THEORY 

In the framework of QCD (including photons), the energy levels and decay widths of 
pionium are functions of the fine-structure constant a ~ 1/137, of the quark masses and 
of the renormalization group invariant scale of QCD. In the following, we concentrate on 
the width F of the ground state. It can be expanded in powers of a and of the quark mass 
difference md — rriu (up to logarithms) . The leading and next-to- leading order terms in this 
expansion are due to the decay into two neutral pions 0, 

F = F2.o+0(5^), 

r2.o = ^aV(ao - a2)' + 0(6'/') ; p'^ = (M^^ - M'^o - \m'^WY^' , (2.1) 
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where oq and 02 denote the S'-wave tttt scattering lengths with isospin 1 = and 1 = 2, 
respectively. We count a and (m^^ — m„)^ as small parameters of order S. The leading 
term in the decay width is then of order 5"^/^. We describe in the present article in detail 
the evaluation of r27rO up to and including terms of order 5^/^, providing details omitted in 
Refs. 
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A. Non-relativistic Lagrangian 

The method used in [|l^] for describing the decay of loosely bound states is an adaption 



of the procedure proposed by Caswell and Lepage some time ago for describing bound 
states in quantum field theories. In the present case, we need to formulate a non-relativistic 
quantum field theory that describes strong and electromagnetic interactions of pions in the 
very low-energy region. The relevant Lagrangian is a rather voluminous object - indeed, it 
contains an infinite number of terms. Fortunately, in the present case, only a small subset 
of that Lagrangian is finally needed. 

The mathematical problem to be solved may be formulated as follows: Construct a non- 
relativistic Lagrangian Cnr that contains all terms needed to evaluate the decay width F 
up to and including terms of order S^^"^. We relegate the construction of this object to the 
Appendixes and y, because the intermediate steps require lengthy calculations, whereas 
the final answer is amazingly simple. Indeed, as already mentioned in |T0|, the following 
Lagrangian achieves the goal: 



NR 



Co + Cd + Cc + C 



^0= Y. - 

i=±,0 ^ 



A 



C 



D 



E 4 

i=±fl 



A2 



+ 



Cs = ciirlnln+n^ + C2[7r^7rl(7ro)^ + h.c] + C3 {tcqTCoY 
+ c^nl A 7rl(7ro)^ + nlnlno A ttq + h.c]. 



(2.2) 



where {u A v) = uAv + vAu, and where A~^ denotes the inverse of the Laplacian. 

The Lagrangian contains explicitly only the pionic degrees of freedom - the sole remnant 
of the photons is contained in the Coulomb interaction described by Cc- The mass param- 
eters M^i coincide with the physical masses of the charged (M^+) and neutral (M^o) pions. 
The role of the low-energy constants (LECs) ci, . . . , C4 is discussed below. 

B. Green functions at a = 



The fundamental objects in the non-relativistic theory are Green functions of the pion 
fields. They are most straightforwardly evaluated with path integral techniques. For in- 
stance, the propagators of the free fields, associated with Cq, read 



5 



= «(0|r7r,(x)7rl(0)|0), (2.3) 

where the vfj denotes a free field. The ie-contribution is generated by a damping factor 
— e / (i^x 7r|(x)7rj(a;) in the action. To ease notation, we always omit this term in the 
following. As is seen from the integral representation ( |2.3|) , the propagator Gatrj vanishes 
for negative times, from where we conclude that the free fields 7f j annihilate the vacuum. As 
a result of this, the Lagrangian Cnr conserves the number of pions. This fact is, of course, 
built in - a term hke e.g. (7ro)^7rQ + h.c. would violate this rule. 

We now discuss Green functions in the presence of interactions, and start the discussion 
for the case where the Coulomb term is absent, a = 0. Again, the relevant Green functions 
may be evaluated in the standard manner through the path integral. First, we note that all 
tadpole diagrams vanish in (split) dimensional regularization, and we adhere in the following 
to this convention. The only corrections to the two-point function are mass insertions, 
generated hj Cd. Summing these up, we obtain 



GnrA^) = (27r)-^ / Jl^^ ^ ; ^,(p) = ^M2. + p2 , (2.4) 



with 

{idt - ^Ml^ - A)GnrAx) = • (2.5) 

Next we consider the four-point functions, relevant for elastic vrvr scattering. To be 
specific, we consider the process 

vr+(pi)7r-(p2) ^ vr+(p3)7r-(p4) • (2.6) 
The corresponding connected Green function is 

Gnr{Pz^Pa]Pi,P2) = j d^xi . . . d^x^ 6-^(^1^1+^2^2-^3x3-^4x4) 

X (0|T7r+(a;3)7r_(x4)4(xi)7ri(a;2)|0),. (2.7) 

Some of the diagrams generated by the interactions are displayed in figure There are 
two classes of diagrams: Mass insertions generated by Cd, and bubbles generated by £5 . 
The perturbative calculation is simply performed by an expansion in the number of loops 
and mass insertions. The reason why this expansion is meaningful is the following. In the 
CM frame = Pi +P2 = (-^"5 0); the elementary "building blocks" to calculate a diagram 
with any number of bubbles are given by the loop integral 

MP") ' 



= !^ (M,.(P° - 2il/40)'^' at D ^ 4 , > 2M,. . (2.8) 

The function Jj is analytic in the complex P° plane, cut along the real axis for P'^ > 2M^i . 
As shown below, the contribution to the scattering matrix element is obtained by putting 
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= 2w+{p), where p denotes the pion three momentum in the CM frame. The loop 
integral is then purely imaginary. In the case where charged pions are running in the loop, 
the integral is of order Ip]^^"^ near threshold. For neutral pions in the loop, it is proportional 
to (M^+ - M^o)^/2 g^i- |-}^g threshold P° = 2M^+. In case that some of the vertices contain 
derivatives (denoted by the full circle in figure [l|c), and/or when mass insertions occur in 
internal lines, additional factors |p| and/or (M^+ — M^^o) appear. As a result of this, the 
expansion in the number of loops and mass insertions is at the same time an expansion in 
IpI and in the isospin breaking parameter — M^o. We conclude that, to calculate the 
scattering amplitude at a given order in the momenta or in the isospin breaking parameters, 
only a finite number of diagrams need to be considered. 

We now to discuss mass insertions on the external lines. These have to be summed up 
in order to generate the correct pole positions at p° = \J + p?. On the other hand, the 
insertions in the internal lines can be treated perturbatively. For a detailed discussion of 
this issue we refer the interested reader to Ref. [Q. The Green function is then of the form 

Gnr{P3,P4;pi,P2) = Y[{uj+{pi) -p^y^R^''^{p3,P4;PuP2)- (2.9) 

i 

The scattering amplitude is obtained from R^'^ by putting all momenta on their mass shell, 

\pO=W+{pi) 



i^'^'^^lpo ( ) = (7r+(p3)7r (p4)out|7r+(pi)7r (P2) in), 



= ^(27r)^5^(pi + P2 - P3 - p,)T^i{p3, P4; Pi, P2) , (2.10) 

with normalization 

(vr+(Pi)|7r+(p2)) = i2nf5\p^ - p,) . (2.11) 

Note that, since the two-point function has residue equal to one, the wave function renor- 
malization constants are unity as well. 

A formula similar to (|2.9| ), ( p. 101) holds for any 2^2 scattering process 

7i\pi)n'ip2)^n\p,)7r"'ip,). (2.12) 

The corresponding relativistic amplitudes are related to the non relativistic ones through 

^i?"'''''(P3,P4;Pl,P2) = 4:{Wi{pl)Wk{p2)Wl{p3)Wn^{p4)}^^^ Tj^R''{p3,P4,PuP2) ■ (2.13) 

In the following, we denote the total and relative momenta by 

P = Pl+P2, P=^(Pi-P2). (2.14) 

Unless stated otherwise, we consider scattering processes always in the CM frame P=0. 



C. The low-energy constants - matching 

We discuss the role of the low-energy constants q that occur in the effective theory. 
We first consider the equal mass case = M^o = M^r, discard the Coulomb interaction 
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Lci and write the corresponding LECs as Ci, C2, C3 and C4. The matrix elements for the 
scattering processes 



TT TT 


— i> TT TT , 






— > 7T 71 , 71 71 - 




ttV- 







(2.15) 

are obtained from the residue of the relevant four-point function, as we just discussed. Each 
contribution consists of a product of loop functions Jj, including vertices with derivatives 
and/or mass insertions. Near threshold, the loop expansion generates a power series in p, 

Tnr = fo + |p|/i + p72 + 0(|p|=^) , (2.16) 

where T/v^ denotes a generic elastic scattering amplitude. The coefficients fi depend on the 
constants q, on the pion mass and on the scattering angle. The threshold amplitude /o 
receives a contribution from the tree graph alone. By use of the relation ( |2.13| ), we therefore 
find that 

leM^cs = T°°'°° , (2.17) 

where T^'^ stands for the relativistic matrix element, evaluated at threshold in the equal 
mass case, in the absence of electromagnetic interactions. 

We have not yet specified what relativistic theory we are considering - the relations ( ^.17| ) 
are true for any of these. Let us consider QCD, and represent the threshold amplitudes 
through the relevant scattering lengths in the isospin symmetry limit rriu = rrid. We then 
have 

SM^Ci = 47r(2ao + 02) , 
3M^C2 = 47r(a2 - Oq) , 

SM^cs = 27r(ao + 202) , (2.18) 

where M^- denotes the pion mass in QCD at niu = md- These relations are true to all orders 
in the chiral expansion. 



D. Matching with the chiral expansion 

There is a second possibility to perform the matching. Namely, one may arrange the 
couplings Ci such that C^r reproduces the chiral expansion of the relativistic amplitude to a 
given order in chiral perturbation theory. To arrive at the relevant expression, it is sufficient 
to work out the chiral expansion of the threshold amplitudes at a given order in the chiral 
expansion and to compare the result with ( |2.17| ). At order p^, the chiral amplitudes are 
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T{n^n ^ TT^TT ) = — , 

T(7r^7r ^nn) = ^ 

T^/^o n Ox s + t + u-3M^ 

i (tT TT TT TT ) = — , (2.19) 



where 



M^ = {m^ + ma)B , B = —\{0\uu\0)\ , (2.20) 

and where F is the pion decay constant in the chiral hmit m„ = = 0. In the isospin 
symmetry hmit a = 0, rriu = rud, the parameter M is further related to the pion mass 
through 

Ml = M^ + O(p^) . (2.21) 

The symbols s, t, u denote the standard Mandelstam variables. It follows that 

1 _ 3 _ 1 , , 

ci = ^ + ---' C2 = -^ + ---, C3 = ^ + ---, (2.22) 

where the ellipses denote higher-order terms in the quark mass expansion. With these values 
of the LECs, the tree graphs of Cis[r reproduce the leading order in the chiral expansion 
of the threshold amplitudes. Similarly, C4 can be related to the momentum dependence of 
T(7r+7r- ^ ttOttO), 



E. Including the Coulomb interaction 

We now consider Green functions at order a, and relax the equal mass condition for 
the pions. There are two classes of diagrams: The first one contains the same diagrams 
as Tattj, but now evaluated at 7^ Mj^o, and with couplings q that depend on a and 
iriu — md, see below. The second class contains diagrams with one virtual Coulomb photon. 
Feynman graphs where the Coulomb photon is attached in such a manner that pions must 
propagate in time in order to connect the two vertices - the self-energy graph is an example - 
all vanish. This is because one may close the contour of integration over the zero-component 
of the photon momentum in a half-plane where there is no singularity in the propagators. 
Since the self-energy diagrams vanish, the mass parameters and M^o in the Lagrangian 
may be identified with the physical masses. The two-point functions for the charged and 
neutral pion field are therefore still given by the expression ( |2.4| ). We now consider virtual 
Coulomb-diagrams that are built from diagrams displayed in Fig. |[ The crosses in the 
figure denote mass insertions. We evaluate the contributions from Figs. ^,c and start 
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the discussion with the Coulomb vertex diagram Fig. ^d, with no mass insertion. After 
integration over the zero component of the loop momentum, the integral to be evaluated is 



The contribution to the scattering amplitudes is obtained by evaluating this expression at 
P° = 2w_|_(p). The result is 

K(P, 2w;+(p)) = -^l^pll. - ^ae, + 0(|p|, - 3) , (2.25) 



where 



4|p| 



is the infrared-divergent Coulomb phase 31 



Next, we consider the two-loop diagram Fig. ^jc, omitting mass insertions. Again in- 
tegrating over the zero components of the loop momenta, the corresponding amplitude is 
expressed in terms of 



2 



d\ 1 d% 



(27r)2'^ J po- 2M^+ - 1iVM,+ |1i - - 2M^+ - li/M,4 

Evaluating this expression at P° = 2w+(p), we find 

Pe(2«;+(p)) = |a(^) + 21n ^ - 1 - m| + 0(|pU - 3) , 

m = /^'^'-'^ - ln47r - r(l)| . (2.28) 

The ultraviolet divergences in diagrams that contain Be are removed in the standard manner 
by adding counterterms to the Lagrangian Cnr- For the consistency of the method it is 
important to notice that the diagrams obtained by adding mass insertions and/or using 
vertices with derivative couplings, are suppressed by powers of momenta with respect to the 
leading terms Be and Vc- They will not be needed in the following. 

The structure of the elastic scattering amphtudes up to and including terms of order a is 
now as follows. First we note that, since the propagators are not affected by the self-energy 
diagrams, the reduction formulae ( p.9| ), ( p.lO| ) are still valid. We write the generic scattering 
amplitude as 

Tnr = T^R + «T^K + 0(«') , (2.29) 

where T^r contains one virtual Coulomb photon. The expansion of the first term in powers 
of the center of mass momenta is as in Eq. ( |2.16| ), with coefficients fi that now also depend 
on the pion mass difference, and on a through the coupling constants q. Omitting the tree 
contribution from one-Coulomb exchange displayed in Fig. we write the second term as 
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^1 

-'-NR 



iPl 



90 



\p\gi + IpI In 



2|p| 



■92 



P'^93 



(2.30) 



The coefficients Qi contain in general infrared divergences, generated by the vertex diagram 
Vc- Otherwise, the structure of the Qi is again the same as the one of the coefficients fi 
in T^fl. Power counting also works in this general case: there is only a finite number of 
diagrams that contribute to a given coefficient fi or gi. Finally, the relation to the amplitude 
in the underlying relativistic theory is again given by Eq. ( p.l3| )p]. 

One may perform the matching to the chiral expansion also in this general case. First, 
we note that a is then counted as a quantity of order p"^. Second, the chiral representation 
( p.l9| ) is valid at order also in the presence of electromagnetic interactions, provided that 
one i) identifies the quantity M in (|2.19| ) with the neutral pion mass, and ii) adds the one- 
photon exchange amplitude in the tt+tt" — > tt+tt" channel. Let us match the amplitudes at 
order p^. Counting powers of F^, it is easy to see that loop diagrams T^^ do not contribute 
-the matching relations become [1^ 



1 



Cl 



C3 



2F2 
1 

16F2 



[1 + k) + 



+ ■ 



C2 



C4 



K 

8F2 ^ 6 



(2.31) 



with K, = M^+/M^o — 1. The ellipses stand either for terms at o(k), or higher-order contri- 
butions in the chiral expansion. The terms of order k are proportional to a at this order in 
the chiral expansion - Eq. (|2.31| ) displays the a-dependence of the couplings q mentioned 
above. 

This concludes our discussion of the evaluation of Green functions in the non-relativistic 
theory. 



III. PIONIUM IN THE NON-RELATIVISTIC FRAMEWORK 

The bound states and their decays are most conveniently described in a Hamiltonian 
framework. The effective theory discussed above renders the pertinent calculations rather 
straightforward, as we will now show. 

A. Hamiltonian and Fock space 

The non-relativistic Lagrangian C^ji gives rise to the following Hamiltonian, 



Both the relativistic and non-relativistic amplitudes must be evaluated by using the same infrared 
regulator, such that the Coulomb-phase can be identified on both sides. We find it convenient to 
use dimensional regularization. 



11 



H = Ho + + He + H5 = Ho + He + V, 
Hr = y"c/^x7^r(0,x); T = 0,D,C,S, 



T^c = 47ra(7rl7r_)A"^(7r|7r+) , 

Hs = -Ci7r^7rL7r+7r_ - C2[7r^7rL(7ro)^ + h.c] - Cg (tt^tto)^ 

- C4[7r:^ A 7r^7ro)2 + ttIttLtto A ttq + h.c] . (3.1) 
It is convenient to introduce creation and annihilation operators, 
[a,(p), al(q)] = (27r)353(p _ q)^^, ■ t,k = ±,0, 

7r,(0, x) = I du{p)e'^-a,ip) ; du{p) = ^ . (3.2) 

The free Hamiltonian becomes 

Ho = / rfz/(p) 5:(m.. + ^)al(p)a,(p) , (3.3) 

and the propagator, evaluated with the free fields 

7ii{x°, x) = e*"«^'°7ri(0, x)e-^"«^° , (3.4) 

of course agrees with the expression Eq. (|2.3|) . We will also need the two-particle states 
with zero total charge, 

|P, p)+ = aUpi) al(p2) |0) , |P, p)o = aS(pi) aS(p2) |0) . (3.5) 

In terms of these, the unperturbed pionium ground state is given by 

1^0, P) = y'dz/(q)*o(q)|P,q)+, (3.6) 

where ^I/o(ct) is the Coulomb wave function in the momentum space 

^o(q) = j^,^, 7 = -^aM^^, (3.7) 

and 

(H„ + He) I "to, P> = (Bo + Bcj;) I *o, P> , 

E, = 2M..{l-^), EcM = ^. (3.8) 

The perturbation V renders the ground state unstable. We discuss in the remaining part of 
this article how the corresponding width can be evaluated. 

12 



B. Resolvents - the master equation 

To determine the width of the ground state, we have considered in Ref. [|T^ the scat- 
tering amphtude in the neutral channel, tt^tt" — > tt^tt^, and determined the position of its 
poles in the complex energy plane. Here, we instead make use of resolvents. While the 
two descriptions are perfectly equivalent, we find that the use of the resolvent renders the 
calculations even simpler. We begin the discussion with the quantity 

GcW = -— (3.9) 

z — tlQ — tic 



whose matrix elements between the charged states ( p. 51) develop poles at the position of the 
energy levels of the unperturbed pionium. To remove the CM momentum of the matrix 
elements, we introduce the notation 

(q|r(2;)|p) = J rfz/(P) +(P, q|R(z) |0, p)+ , (3.10) 

where R(-z) denotes any operator in Fock space. One can now easily relate the matrix 



element of to Schwinger's Green function |32 



(27r)353(q_p) 1 4™ 1 

(q|gc7(2:)|p) - 



E - qVM,+ E - qVM^+ |q - E - pVM^+ 

^ A7rar]I{E; q, p) \—— , (3.11) 



E - qVM,+ ' ' ' ^' E - pVM, 
with 

/(E;q,p) = y^ [(q _ ^ r]ya^{l - xY{E - qyM^+){E - p7M,+ )] ' ^^"^^^ 

where r] = {—E /M^+)~^^'^ and E = z — 2M^+. This function has poles at = 1, 2, 

In order to calculate the position of the poles in the real world, with V 7^ 0, we consider the 
full resolvent 

G{z) = ^—. (3.13) 
z — ti 

Expanding in powers of the perturbation V, one finds that G satisfies the equation 

G = Gc + Gc-tGc , 

T = V + VGcT. (3.14) 
We remove the ground state singularity from Gc, 

Gc = Gc{l- y'MP)l^o,P)(^o,P|} , (3.15) 

introduce 
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r = V + YGcT , 



(3.16) 



and find for G the representation 

G = Gc + GctGc + (1 + GcT)no(l + tGc) , (3.17) 



where 



EcM - Eo - {'^o\f{z;P)\^ 



and 



(^o|t(z;P)|^o) = / MP')(^o,P'|r(z)|^o,P). (3.19) 

The singularity generated by the ground state pole is absent in the barred quantities. There- 
fore, the pertinent pole must occur through a zero in the denominator of the expression 
( p.l8| ). In the CM frame, the relevant eigenvalue equation to be solved is 

z-Eo-(^o|t(z)|^o) = 0, (3.20) 

where the matrix element denotes the quantity on the left-hand side of Eq. ( |3.19|) , evaluated 
in the CM frame P = 0. 

The master equation (|3.2CI|) is a compact form of the conventional Rayleigh-Schrodinger 
perturbation theory. Note that it fixes the convergence domain of the perturbation theory: 
the theory is applicable as long as the energy-level shift does not become comparable to 
the distance between the ground-state and the first radial-excited Coulomb poles. Equation 
( p.20| ) is valid for a general potential- containing e.g. the interaction with the transverse pho- 
tons - since in the derivation, we did not use the explicit form of the interaction Hamiltonian 
in (K 



C. Singularity structure of the resolvent 

We find it instructive to shortly discuss the analytic structure of the matrix elements of 
the resolvent G, and the location of the shifted ground state pole. First, from ( |3.20| ), it is 
seen that this pole will occur at the same position for any channel. Second, it is expected 
on general grounds that the pole will move to the second Riemann sheet. Indeed, consider 
the operator r in the second iterative approximation 

r = V + VGcV + 0(V^). (3.21) 

To evaluate the matrix element between charged states as required, we insert a complete set 
of neutral states in the second term. The eigenvalue equation becomes 

z = Eo + {-ci - 2clJo{z) + ■••}, (3.22) 

where Jq denotes the loop integral (|2.8|) . This function has a branch point ai z = 2M^o, and 
its imaginary part has the same sign as the imaginary part of z throughout the cut z-plane. 
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Therefore, the equation ( p.22|) has no solution on the first Riemann sheet. On the other 



hand, if we analytically continue Jq from the upper rim of the cut to the second Riemann 
sheet, we find that a zero at 



z = Kez + i Imz 



Rez = Eo-^^Ci + --- , lmz = T^^P^ + - , 3.23 

with p = 2M^o(M^+ - M^o - 

The imaginary part is of order 5''/^. We demonstrate below that it is the only term at 
this order. Using Eq. ( |2.18| ) and F = — 2Im2;, one recovers the leading order result 



Similar arguments apply to all the other pole positions. [Of course, in order to correctly 
describe the new positions of the exited energy levels, our original Lagrangian Cnr must be 
enlarged.] We conclude that the 2-particle matrix elements of G are analytic functions in 
the complex z-plane, cut along the real axis for Re^; > 2Mt^o. The poles are located on the 
second Riemann sheet. 



IV. PIONIUM DECAYS 



A. Perturbative solution of the bound-state equation 



In order to find the solution to Eq. (|3.20|) at order (5^/^, it is convenient to reduce the 
equation ( ^.16| ) to a one-channel problem with an effective potential W. We use a projector 
g on the two-particle states |P, q)+. 



f/z/(p)rfKq)|P,q)++(P,q| 
1 - ^, 



(4.1) 



and find in the standard manner 



QTQ = gWg + gWg Gc Qtq , 

W = V + \q,Gc {l - f?oVf?oGc}"' QoV . (4.2) 

This is result is still perfectly general. In the case considered here, one may simplify the 
expression for the effective potential, replacing Gc by Gq = {z — Ho)~^ 

W = V + Vf?oGo {1 - QoVqoGoY' f?oV . (4.3) 

The matrix elements of the effective potential can be expanded in powers of momenta, 
because there are no nearby singularities. Specifically, we write 

(q|w(z)|p) = {27^f5\c^ - p) (--^ + ■ ■ ■) 

7r+ 

+ w{z) + Wi{z)p^ + W2{z)c^ + W3(z)pq H . (4.4) 
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If we now iterate the equation ([4.2|), at the order of accuracy we are working, the decay 
width of the 7r~^n~ atom r27rO = — 2Im z is given by 

T2n0 = ^ lmw{l + 2Rew (gc(^o)) ) + 0{6') , (4.5) 

where w = w{Eo), and 



{2ttY (27r)' 

In order to calculate this integral, one needs to define Schwinger's Green function in d 
dimensions. In field theory, the Fourier transform of the Coulomb potential in d dimensions 
is given by exactly the same expression as in 3 dimensions - consequently, the first two 
terms in the representation ( p.ll| ) are also valid at d 7^ 3. For the last term, the integral is 
convergent, and we may work at d = 3. The integral is then equal to 

{gc{Eo)) = ^ e , e = 2 In a - 3 + A(/i) + In ^ . (4.7) 



The quantities Re w and Im w can be determined from iterations of Eq. ( [4. 21 ) to the needed 
accuracy. 

Re w = — ci , 

to„ = -^P"^ (1 + - 2pc,)^(l - p^) . (4.8) 

We have now expressed the decay width of the tt+tt^ atom in terms of the non-relativistic 
couplings Ci ■ ■ ■ C4. It remains to determine the relevant combination of these couplings from 
the matching of the relativistic and non-relativistic amplitudes. 



B. Matching to the threshold ampHtude 

We determine the non-relativistic couplings that enter the expression for the decay width 
through Eq. ([4.8|) , and start the discussion with the couplings C\ and C3. These contribute 
to the decay width at order 5^/^, because p counts as a quantity of order b. Therefore, these 
two couplings are needed at order 5" (no isospin breaking), as a result of which we may 
replace them by the isospin symmetric quantities c\ and C3 in Eqs. ( |2.18| ). It remains to 
determine the combination of the couplings C2 and C4 that enter ([4.8|) . As we will now show, 
it suffices for this purpose to calculate the real part of the tt+tt" — ^ 71*^7?° amplitude at order 
5 in the non-relativistic and in the relativistic theories. 

Whereas the Hamiltonian framework is very convenient to discuss the energy spectrum, it 
is more convenient to calculate scattering amplitudes in the Lagrangian framework discussed 
in section ||. In the relativistic theory, the on-shell amplitude for vr+vr" tt^tt^ contains 
infrared singularities that exponentiate P5 |, 
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^00;± _ ^^^(^, dOO;±N mOO;± oOO;± _ ( ^^(0 



where pi, p2 denote the 4- momenta of incoming vr^ and vr^ mesons. In Ref. I^Sj it is 
demonstrated that - using a photon mass as an infrared regulator - the residual amplitude 
T^'"^ is free of infrared singularities. Here we assume that the same is true in dimensional 
regularization. We find 

= z^, + ^ - f + 0(|p|, D - 4) , (4.10) 

4|p| /TT 



where 9c is defined by Eq. ( p.26| ). The infrared divergences cancel in the real part of 5°°'='= 



at threshold, whereas the imaginary part is divergent at Z) — 4. 

One may verify that at order a, exactly the same divergent Coulomb phase appears in 
the non-relativistic amplitudes. Indeed, if one performs the calculation at 7^ 3 and splits 
off the phase according to 



T;;r(q,p) = e^"'^^T;;^"(q,p), (4.11) 

then there are no infrared singularities in the amplitude T^'^'^(q, p) at threshold in the limit 
(i — s> 3, at order 5. For the real part, we findQ 



Ref'^i{^. P) = § + 1^ Ir^ + I^iT^ ^^^t1.7°° + «(P) , (4.12) 

where 



Bi = C2 + o{5) , B2 = C1C2 + o{5) . (4. 13) 

The singular contributions ~ l/|p|,ln|p| are generated by the exchange of one Coulomb 
photon (see Figs. Hb,|^). At 0(5), the constant term in Eq. (|4.12|) is equal to 

1 „ „ ...2 ( c^dM'n\ . cM' /■ A.f 



Re = 2c, - 4Af,.n{c, + + ^ {l - A(a) - !« ^) c,c„ (4.14) 



where A(/i) is given by Eq. (|2.28| ), and k, is the same as in Eq. ( |2.31|) . The ultraviolet 



divergence contained in A(/x) may be absorbed in the renormalization of the coupling C2- 
This procedure at the same time eliminates the ultraviolet divergence in the expression for 
the decay width. 

In the following, we assume that - up to and including terms of order S - the relativistic 
amplitude does have the same singularity structure - as a function of the momentum p - 
as the non-relativistic amplitude ( ^.12[ ). We can then match the non-relativistic expression 



Note that in Ref. @, the non-relativistic scattering amplitude is defined with an opposite sign. 
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0° in 



to the relativistic one in the standard manner, using Eq. ( p.l3| ). The quantity ReAj^^ 
( [4.12| ) corresponds to the one introduced in Ref . |^| , in the context of the relativistic theory 
(modulo the Coulomb phase, which does, however, not contribute to the amplitude at order 
The logarithmic singularity is absent in the amplitude at order e^p^ at which the 
were carried out - it first emerges at order e^p^ , see Appendix 0. 



calculations in Ref. 

Finally, the relation ( [4.14| ) represents the matching condition between the regular part of the 
relativistic vr+vr" tt^tt" scattering amplitude at threshold and the pertinent combination 
of non-relativistic coupling constants q. 



C. General expression for the decay width 

Substituting the results of the matching into the expression for the tt+tt^ atom decay 



width ([4.5[), and using (|4.7|) and (^78|), we obtain 



= ^ aYA\l + K), ^ = - J- Re A+7°° + o{6) , 

K = - {ao + 2a2f - — (Ina - l)(2ao + 02) + o{S) , 
y o 

= (M^ - M^^o - iM'^.aY' ■ (4.15) 

This is the general expression for the tt+tt^ atom decay width, valid at next-to-leading order 
in isospin breaking, and to all orders in the chiral expansion. Note that all mention of the 
non-relativistic theory has disappeared in the final result that relates the observable quantity 
(the decay width) to the relativistic scattering amplitude at threshold. 

The primary objective of the DIRAC experiment is to measure the difference ao — a2 
of the S-wave tttt scattering lengths that are defined in the isospin-symmetric world. The 
expression ( |4.15|) is not yet suited for this purpose, because it relates the width to the 
scattering amplitude at threshold. This quantity contains the combination oq — 0,2 we are 
looking for, together with isospin breaking contributions. One has to evaluate these and 
subtract them from the measured amplitude. ChPT allows one to achieve this goal - order 
by order in the expansion in the quark mass. 



D. Amplitude at O(e^p^) 

The normalization of the quantity A is chosen such that, in the isospin symmetry limit, 
it coincides with the difference — ^2 of the S-wave scattering lengths. In the general case, 
we expand the amplitude in powers of the isospin breaking parameters a and — rrid, 

A = aQ — a2 + hiim^ — rriuY + h2a + o(5) . (4.16) 

This decomposition is true irrespective of the chiral expansion. The scattering lengths as 
well as the coefficients hi are functions of the quark mass m and of the renormalization group 
invariant scale of QCD. What is the meaning of ag — 02 in the presence of isospin-violating 
interactions? To clarify the issue, we consider the expression A at leading order in the chiral 
expansion. From (|2.19|) , we find 
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A = ^^(4M.'+ - M') + 0{p\ eV) . (4.17) 

To bring this into the form ( [4.16| ), we note that, in the isospin-symmetry hmit m„ = rud, a = 
0, the scattering lengths can be expanded in powers of the pion mass, defined to be the 
position of the pole in the correlator of two axial currents. It is an algebraically perfectly 
legitimate procedure to identify this mass with the charged pion mass. We adhere in the 
following to this procedure, in order to agree with the standard conventions in ChPT. The 
expression for the difference of the scattering lengths then reads 

ao-a2 = ^^ + 0(p^). (4.18) 



Comparing this with ( |4.17| ), we find 



A = ao-a, + + 0{p\ eV) . (4.19) 

From this result, we may read off the coefficient /i2 at leading order in the chiral expansion, 

32a7rF2 ^ ^ ^ ^ 

[To be precise, the first term on the right-hand side of this equation should be evaluated at 
a = 0. To ease notation, we omit this request here and in the following]. On the other hand, 
the above calculation is not accurate enough to determine hi at leading order, because for 
this purpose, the amphtude is needed at order p^. This procedure may obviously be carried 
out order by order in the chiral expansion - all that is needed is the chiral expansion of the 
scattering amphtude at threshold, at niu 7^ rria, a ^ 0. As a result of this, the quantities hi 
are represented as a power series in the quark mass m (up to logarithms). 

The evaluation of the amplitude for vr+vr" —* tc^tt^ has been carried out at 0{p'^,e^p'^) 
in Ref. This result allows us therefore to determine the coefficient hi (/i2) at order j9° 



(p^). Some remarks are in order. 

i) In Ref. [^], the scattering amplitude has been evaluated at m„ = m^. For our pur- 
poses, the expression for generic and nid is needed. On the other hand, up to and 
including terms of order p^, the strong amplitude does not contain mu — md terms. The 
only source for such contributions is the tree graph, where is expressed in terms of 
the neutral pion mass. The generalization of the result Ref. to the unequal mass 
case is therefore straightforward. 



ii) The normalization point in Ref. |3J] is chosen to be the neutral pion mass. According 
to our definitions, we have to normalize all low-energy constants at the charged pion 
mass. The terms that emerge from the shift of the normalization point are proportional 
to = — M^o and are included in the expressions given below. 

The rest is then straightforward. We find 



hi = 0{m) 

h2 = IZTT-^ I 1 + 



23 - 3 - 
hli + - k 



3 M 2 

-pih) +0(m2), (4.21) 



2567r2F2 
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where p{ki) stands for the following combination of the electromagnetic low-energy con- 
stants [H], 



p{ki) = -30 + 9ki + 6k3 + 2ke + k8 + - Z{ki + 2k2 + Gk^ + Uka - 6ks) 



(4.22) 



and 



A! 



A. 



A em 

TT 



(4.23) 



The quantities ki denote again the running coupling constants A;[(/i) at scale fi = Mt,+ . Note 
that according to our counting, the quantity = ln(M^+/M^o) introduced in Ref. p4 



IS 



of order S and hence does not contribute to h2- Further, F may be expressed through 
according to M 



F 1 + 



Ml 



167r2F2 



(4.24) 



For the numerical analysis, one has to specify the values of the low-energy constants that 
enter the expression for h2- We are not aware of an estimate for the SU{2) x SU (2) couplings 
ki. On the other hand, the corresponding couplings Ki in SU{3) x SU{3) |^ have been 
estimated by invoking e.g. sum rules or a resonance saturation hypothesis ||37|-p9|. In order 
to use this information, we need to relate the couplings ki to their SU{3) counterparts Ki. 
In Appendix |^, we show that 



where 



1287r2 



p{ki) = P{Ki) - 8Zh 



-6(K[ + K;) + SKI - 5K; + KI + Q{Kl + KI, + KI,) 



(4.25) 



;i8 + 28Z)ln^-2zfln^ 



30 



(4.26) 



and where K^ denote the running couplings introduced in Taking into account this 

relation, we may rewrite the formula for the width in the following form. 



with 



12^0 = -aV (ao - 02 + e)\l + K) 

y 



127r2F2 



23 - 3/ 

y + ^^ + 4^^. 



32vrF2 

-|-0(m(m„ - m^)^, arr?) o(5) . 
The quantity K is given in Eq. (|4l5| ). 



3aM2+ 
^ 256^^ ' 



(4.27) 



(4.28) 
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E. Numerical analysis 



In the numerical evaluation of the lifetime, we use for oq and 02 the values from the 
recent analysis in Ref. 0, Oq = 0.220 ±0.005, 02 = -0.0444 ±0.001, 09-02 = 0.265 ±0.004. 
To evaluate the correction e, we first recall that the non-electromagnetic part of the pion 
mass difference is tiny, of order ~ 0.1 MeV [^. Therefore, we identify A^™" with the 
experimentally measured total shift A^. Further, in the calculations we replace ttIsBq by 
M|^+ — M^+/2, according to our definition of the isospin symmetry limit. The values used for 
the low-energy constants in the strong sector are li = — 0.4±0.6, Z3 = 2.9±2.4 0]. For KKn), 
we use the values given by Baur and Urech in Ref. [^, Table 1]: K[ = —6.4, = 6.4, = 



-6.2, = 19.9, = 8.6, = K{q = 0, = 0.6 (in units of lO'^). We evaluate 
P{Ki) at scale n = Mp. Further, we attribute an uncertainty 2/167r^ - that stems from 



dimensional arguments - to each Kl. The values of Kl obtained both by Moussallam |38 



and by Bijnens and Prades lie then within the uncertainties attributed. The same is 
true, if the saturation is assumed not at scale /i = Mp, but somewhere within the interval 
0.5 GeV < fi < 1.0 GeV. Finally, we use = 92.4 MeV. Adding the uncertainties in 
Zi, Z3, ao, a2 and in Ki quadratically, we obtain 

e = (0.61 ±0.16) X 10"^ = (1.15 ±0.03) X 10"^ (4.29) 



or 



2 

9 



2^0 = -aV (00-02)^(1 ±5r), 5r = (5.8 ±1.2) X 10-2. (4.30) 



This amounts to a six percent correction to leading-order formula by Deser et al [0]. In 
the total decay width F, the decay into 27r° is by far the dominating mode. For example, 
the decay width into a 27 pair, which is the first subleading mode in 5 counting, is T2-f = 
Mt,+ /A [|l],0 at leading order in the 5-expansion, as a result of which one has r2^/r27r0 ~ 
3 X 10^'^. For this reason, one may safely identify T2t^o with the total lifetime, 

r = F-Jo = (2.9±0.1) X 10"^^ s. (4.31) 

We add the following remarks concerning these numbers. 

• The bulk part in the uncertainty in the lifetime is due to the uncertainty in the 
difference of the scattering lengths Oq — (12-, which results in ±0.085 x 10"^^ s. 

• The uncertainties in the constants increase this to ±0.091 X 10"^^ s. Including the 
remaining uncertainties does not change this number in the digits displayed. 



The numbers in Eqs. (|4.29|) - (|4.31| ) differ from the corresponding ones in our previous 



paper [12], because the present values of the scattering lengths, of li and of differ 



from the ones used there. The above values of ao,a2 are the result of a complete 



analysis at order p - they replace the ones used in |12], taken from the preliminary 



numerical result cited in Ref. |42|. The present value of Zi is based the same analysis [Q. 
The bulk part in the change of the lifetime is of course due to the updated value of 
00 — 02, because this combination of scattering lengths enters the expression for the 
decay rate at leading order. 
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• The vacuum polarization correction to the hfetime, that is not taken into account here, 
amounts [^] to a contribution of —0.01 x 10~^^ s. 

We expect that the higher-order contributions to the vt+tt^ atom decay width in ChPT 
are neghgibly small. Consequently, an accurate determination of oq — 02 from a precise 
lifetime measurement is indeed feasible. 



V. SUMMARY 

i) We have considered decays of the tt+tt^ atom in its ground state. Aside from a 
kinematical factor, the decay rate can be expanded in powers of the isospin breaking 
parameters a and (m^ — m^^)^. It is convenient to book these parameters as terms of 
order b. 

ii) To calculate the leading and next-to-leading order term in the (5-expansion of the 
width, we have constructed a non-relativistic Lagrangian that describes the low-energy 
interactions of pions and photons. In this framework, the matrix elements of the 
resolvent 1/(2 — H) develop poles on the second Riemann sheet in the complex z-plane. 
The positions of the poles are related to the energy levels and widths in the standard 



manner. By using Feshbach's technique, we have derived the master equation (|3.20|) 
for the position of the ground-state pole. 

iii) On the basis of this equation, we have calculated the decay width of the ground state 
of pionium in terms of the parameters of the non-relativistic Lagrangian. At leading 
and next-to-leading order in the 5-expansion, only the channel A^+^- tt'^tt'^ is open. 
Furthermore, at this order of accuracy, transverse photons do not contribute - the 
relevant Lagrangian becomes then very simple, see Eq. ( ^.21) . Matching the non- 
relativistic amplitude to the relativistic one, we have then expressed the decay width 
in terms of the relativistic scattering amplitude, up to terms that vanish faster than 
(5^/^. The relevant formula is displayed in Eq. (|4.15| ). 

iv) At this stage, one may invoke ChPT, which allows one to expand the isospin breaking 
terms in powers of the quark mass, and thus to get contact with measurable quantities. 
The result is given in Eq. ( [4.27| ), that displays the width in terms of the combination 



ao — (32 of S'-wave scattering lengths, and a correction that we have calculated at order 
a and ma. The quark mass difference shows up only at order (m„ — md)^m. We 
expect this term to be completely negligible. The recently determined values of the 
TTTT S'-wave scattering lengths gives 

r = (2.9±0.1) X 10-^^s. (5.1) 

Since the isospin breaking corrections at order a and ma are small, we expect that 
chiral corrections at higher order as well as higher-order terms in isospin breaking are 
irrelevant for data on the lifetime obtained in the foreseeable future. 
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APPENDIX A: GENERAL NON-RELATIVISTIC LAGRANGIAN 

In this Appendix, we outline general rules for the construction of a non-relativistic La- 
grangian that describes low-energy interactions of pions and photons. The Lagrangian does 
not contain terms that correspond to transitions between sectors with different number of 
heavy particles (pions), since these belong to hard processes and are hidden in the couplings 
of the non-relativistic Lagrangian. For this reason, in order to describe tttt scattering in 
the non-rclativistic framework, it suffices to consider Lagrangians in the sectors with one or 
two pions (including any number of photons). The theory must be invariant under space 
rotations, C, P, T and gauge transformations. On the other hand, due e.g. to the presence 
of photons, the Lagrangian is not invariant under Galilei transformations. The appropriate 
building blocks are provided by the covariant derivatives of the charged pion fields 

DtT:± — dtTT±^ ieAQT:± , D 7r± = V 7r± ± ie A 7r± , (Al) 

and the electric and magnetic fields 

E = -VAq - A , B = rot A . (A2) 

For the neutral pion field, the covariant derivative coincides with the ordinary one. 

The Lagrangian consists of an infinite tower of operators with increasing mass dimension. 
All possible operators allowed by the symmetries must be included. In particular, in the 
one-pion sector, the Lagrangian is given by 

£. = i(E^-BV4{>9.-M,. + ^ + ^ + ...}.„ 

+ i;-i{«.-M,. + ^ + ^ + . (A3) 

complemented, e.g., with all possible non- minimal couplings containing E and B, that we 
have not explicitly displayed. At tree level, this Lagrangian reproduces the relativistic result 
for the scattering amplitude tt -|- 777.7 — >^ tt -|- 777 at 0(e"^"'""), to all orders in the momentum 
expansion. 

In the two-pion sector of zero total charge, one has to construct the operators that contain 
four pion fields and any number of photon fields. The lowest-order Lagrangians with zero 
and two space derivatives are given by 
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= CitiXti^_ti+ti_ + C2(7r;^7rl7ro7ro + h.c.) + CavrJ^TToTro , (A4) 
= C4|(4 7rl)(7ro7ro) + (4^-)K d' ttq) + h.c.j 

+ C5|(7r^ 7rl)(7r+7r_) + h.c.j + C6|[7r^7r+] B^ [7rl7r_]| 

+ crj (4 B^ 4) (ttotto) + h.c. I , ( A5) 



»2 



where uY)v = u D^f + f D^u. Note that the couphngs are not necessarily real, as a result 
of which the Lagrangian is not, in general, hermitian (see below). Again, we do not display 
explicitly non-minimal couplings that, apart from covariant space derivatives, contain the 
vectors E and B. Moreover, we do not display covariant time derivatives, or higher-order 
Lagrangians £2^'*, c!2\ ''' which contain 4, 6, ■■■ space derivatives. In the absence of 
photons, the Lagrangian £2°"* + ^^2'' is equivalent to the one given in Ref . . 

In the isospin symmetry limit a = 0, = m^, the following relations hold for the 
corresponding couphngs Ci, ■ ■ ■ C7, 

C3 = ^ (ci + C2) , ^7 = ^ (c4 + C5 + ce) . (A6) 

In the sectors with one or two pions, the Lagrangian is therefore given by 

£ = £1 + + £f + • • • , (A7) 

where the ellipses stand for non-minimal, or higher- dimensional operators or higher-order 
terms with covariant time derivatives. On the mass shell, the latter terms are eliminated 
by using the equation of motion (EOM). However, they need to be included if one decides 
to renormalize Green functions At tree level, the Lagrangian ( [A7| ) reproduces the 

relativistic result for the scattering amplitude /ctt + rwy ^ k-K + wj, k = 1,2 at 0{e^~^^), to 
all orders in the momentum expansion. 

The scattering amplitudes in the non-relativistic theory are related - through the reduc- 
tion formula - to the residues of the pertinent Green functions in a standard manner. [Note 
that, in the non-relativistic theory, one has to sum up all insertions ~ p^, p^, ■ ■ ■ in the 
external legs of the pions, in order to ensure that the poles sit at the correct position, ac- 
cording to the relativistic dispersion law p° = (M^j + pf)^^^. On the other hand, insertions 
in the internal lines are treated perturbatively. Fore more details, we refer the reader to 
Ref. [0.] The loop corrections to the Green functions are then calculated in a standard 



manner, by using Feynman diagrammatic techniques, with one important modification. It 
is well known that in the non-relativistic theory, in the presence of light particles (photons), 
the Feynman integrals should be properly butchered, in order to avoid contributions from 
the loop momenta at a hard scale - otherwise, loop corrections to the Green functions would 
lead to a breakdown of the counting rules in the non-relativistic theory. A suitable procedure 
built on top of the Feynman rules in the non-relativistic theory is provided by the so-called 
threshold expansion |Q (see also [Q), that enables one to disentangle the contributions 
coming from different regions of loop momenta, by expanding the integrands - in the dimen- 
sional regularization scheme - in all possible small kinematical variables. Here, we adopt 
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a simple and physically transparent formulation of such a procedure [^]. First, we put a 



momentum cutoff on all 3-dimensional Feynman integrals, after having performed the inte- 
grals over all zero components of virtual four-vectors (by eventually using split dimensional 
regularization). Then we choose the cutoff mass to be much smaller than the hard scale, 
given by the pion mass. Next, we expand the integrands both in external and integration 
momenta. In the presence of the cutoff, this is a well-defined procedure. At the last step, we 
remove the cutoff and calculate the integrals, with the expanded integrands, in dimensional 
regularization. This sequence of steps systematically removes the hard-momentum contri- 
bution from the integrals, which at low energies is given by a polynomial in the external 
momenta. 

All couplings in the non-relativistic Lagrangian are determined by matching to the rel- 
ativistic theory. In the presence of photons, the non-relativistic Lagrangian is, in general, 
not hermitian - because these constants are not real. This is due to the fact that in the 
non-relativistic approach, one has shielded some of the intermediate states - those with 
masses below the two-pion threshold - that appear in the relativistic theory and that belong 
to the class of hard processes in the non-relativistic terminology. The imaginary part of 
such diagrams then contributes to the imaginary part of the couplings of the non-relativistic 
Lagrangian. For example, the decay of the 7r"'"7r~ atom into two photons in the relativistic 
theory is described - at leading order in a and in rn - by the imaginary parts of the diagrams 
depicted in Fig. |^. These diagrams are not present in the non-relativistic theory. On the 
other hand, they contribute to the imaginary part of the coefficient c\ at order 0[e^\ 

The next remark concerns power counting. In fact, we have three different types of power 
counting in our theory: 

i) Non-relativistic power counting. The Green functions - calculated at a fixed order 
in an expansion in the coupling e - are expanded in powers of external 3-momenta 
of the particles, and in the mass difference — M^o in the manner described in 
Sections fTE] and |ITE| . 

ii) Chiral power counting. After matching to ChPT, the couplings are given in a form of 
a series in the quark masses and e. The coefficients contain the low-energy constants 
(LECs) of ChPT. This procedure is systematic in the sense that matching at higher 
chiral order does not affect the result obtained at lower orders. 

iii) Counting the isospin-breaking parameter b. After matching to ChPT, the couplings 
in the non-relativistic Lagrangian can be rewritten as 

a = cf ) + acf) + (m, - m^^cf + o{6) , cj") = (M^^ , LECs) , (A8) 



where, by convention, we have defined the isospin-symmetric world with a = 0, and 
nT'd = fTT-u as the one in which the common mass of the pion triplet coincides with the 
mass of the charged pion in the real world [Note that in the relativistic pion scattering 
amplitudes, odd powers of — niu never appear.]. Consequently, the powers of 6 in 
the expansion of the Green functions around the isospin symmetry limit stem from 
different sources. The explicit powers are due to the coupling to photons, and the 
implicit powers are encoded in the couplings of the Lagrangian, as well as the charged 
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and neutral pion mass difference. In the calculations, one has to carefully keep track 
of all these sources of corrections in a given order in S. 

The Lagrangian { \K^ ) contains an infinite number of operators. In actual calculations, 
only a few of them are needed. In particular, we will make it plausible that, in order to 
calculate the pionium decay width at 0(5^/^), it suffices to work with the Lagrangian given 
in Eq. (|2.2|) . The arguments in favor of this simplification are provided in the following two 
Appendixes. 



APPENDIX B: THRESHOLD EXPANSION AND THE ROLE OF TRANSVERSE 

PHOTONS 

We illustrate the evaluation of the scattering amplitude in the non-relativistic theory in 
the presence of photons. We work in the Coulomb gauge, which allows a clear-cut separation 
of Coulomb and transverse photons, and argue that the radiative corrections to the tt+tt^ — >■ 
7r°7r° scattering amplitude, generated by transverse photons, vanish at threshold at order e^. 



1. Pion self-energy 

We start with the two-point function of the charged pions. The self-energy correction at 
O(e^) due to the diagram Fig. |^a is given by 

.J ..e- (0|T.,(.).t (0)10) ^ ^^,^p./(,^j)_^o_^,o,p) > (Bl) 

' - M^^ J {2nr^ -P{M^, + (p - l)V(2iV4.) - pO + /o) + ^ ^^2) 

The threshold expansion of the above integral amounts to (we place a hat above the 
threshold-expanded quantities) 

^.r. N f d'^l / ^ (pl)2\ 1 r 1 /pi P \ 1 



2M^+J {2nYV P VM^+ 2M^+J {Q+\\\) 

^ = ^-+ + ^-P°' d = D-l. (B3) 
Performing the remaining integration, we obtain 

m p) - p^-.'^-^ r(^)r(2-c^) 



-,2 



p^^] |L(/i) + ln 7Tj + ---, when d^3, (B4) 

m = /"'(^ - I (r'(l) + ln4vr + 1)) . (B5) 
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As usual, /i denotes the scale of dimensional regularization. 

In order to remove the divergence from the two-point function, one introduces the coun- 
terterm 

A/: = E /i4d^(^A - M.. + , /i = L{i^) + rM . (B6) 

The contribution to the two-point function is displayed in Fig. The quantity /[(/x) 
denotes the finite, scale dependent part of the coupling constant /i. 

To calculate the wave function renormalization constant for charged pions, one has to 
reverse the limits Namely, we perform the limit — > (mass-shell limit in the non- 
relativistic theory) at > 3, in order to avoid the infrared singularity. Since the ratio 
E(p°,p)/i7 vanishes as f2 -h> 0, the self-energy diagram Fig. ^ does not contribute to the 
wave function renormalization constant. The sole contribution comes from the counterterm 
given by Eq. (p6D, 

Note that Z±{Q) = 1. This feature is due to the derivative coupling of transverse photons, 
and to the use of the threshold expansion, which guarantees that the non-relativistic power 
counting is not altered by loop corrections. 



2. Scattering amplitude vr+vr —>■ vr'^vr'^ at order 

We now discuss the radiative corrections to the 7t~^tt~ vr'^vr^ scattering amplitude at 
order e^, due to transverse photons. We have to consider diagrams with any number of 
strong bubbles [including of course the tree diagrams], and attach one virtual photon line 
to these diagrams in all possible ways. The photon couples to two pions (the relevant part 
of the Lagrangian is given by Eq. (|A3|)), as well as to the vertices with four pions (see 
Eq. ( |A5| ) ) . Further, the photon couples to two pions in a minimal way, as well as through 
the non-minimal vertices which contain more space derivatives acting on the fields. 

Some preliminary remarks are in order. After applying the threshold expansion to a given 
diagram, one always ends up with a homogeneous integrand, and naive power-counting is 
restored. Since a strong bubble introduces a suppression factor in a diagram with no photons 
(see Section we expect that - even in the presence of photons - diagrams containing more 
strong bubbles will be more suppressed, and, for a given topology, it suffices to consider 
diagrams with a minimal number of strong bubbles. The same consideration applies to 
diagrams with non-minimal photon couplings, and to diagrams with derivative couplings in 
strong four-pion vertices: since power-counting holds, we expect that these are suppressed 
with respect to the diagrams of the same topology, but with a minimal number of derivatives 
in the vertices. 

We start with the diagrams where the photon couples only to two pions. According to 
the above discussion, we do not consider diagrams with non-minimal couplings, and restrict 
ourselves to the non- derivative strong Lagrangian (|A^). The set of all topologically distinct 
diagrams with one virtual transverse photon coupled in a minimal way to two pions, is 
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depicted in Fig. |^. In each class, we single out a representative with a minimal number of 
strong loops. 

The corrections to the external legs (Fig. ^) vanish at threshold, because Z±{0) = 1. 
Next, we consider the diagram corresponding to the exchange of a transverse photon between 
the initial vr^^vr" pair (Fig. Qd). The integral to be calculated in this case is given by 

r d^'l / 2 (pl)2\ 1 



^+-7(IPI) = "XF^ / 77^ (P 



iV47 J (27r)«z V" F y (M.+ + (p-I)7(2M.+ )-pO 

^ (B8) 



/2(M.+ + (p - 1)7(2M^+) - - P) ■ 

One has still to multiply this integral with the purely strong amplitude in order to get the 
contribution of the diagram in Fig. In this expression, p denotes the relative momentum 
of the TT^7c~ pair in the CM frame, and p^, are the energies of tt"*" and tt" particles. 
We put the external particles on the mass shell, p° = p° = + p^/(2M^+) + O(p^), 
and perform the threshold expansion in the integral. Note that with this procedure, the 
integrands also should be expanded in the O(p^) remainder of p^.. The threshold-expanded 
integral in Eq. (p8| ) can be rewritten in the following manner, 

? /, ,x e2 r d'^l 1/2 (pl)'\ 1 

e^pl (B9) 



16M^+ 87rM^+ V ^ 

Again, it is seen that this particular contribution vanishes at threshold. 

We have also investigated the remaining contributions depicted in Fig. ^. All these 
contributions vanish at threshold. Moreover, we have considered all topologically distinct 
sets of diagrams where the virtual photons couple to four-pion vertices, depicted in Fig. 
Again, in each set we have restricted ourselves to the diagram with a minimal number 
of strong loops, and with a minimal number of derivatives in strong and electromagnetic 
vertices. We have found that the contributions from all these diagrams vanish at threshold. 

To conclude, we have considered all topologically distinct diagrams for the scattering 
process tt+tt^ tt^tt*^ in the non-relativistic theory, where one virtual photon couples in 
all possible ways to strong diagrams. From each class of diagrams, we have singled out the 
representative with a minimal number of strong loops, and a minimal number of derivatives 
in the vertices. We have checked that each such diagram vanishes at threshold. Using power- 
counting, the same is seen to be true for the diagrams with more loops, and/or higher-order 
couplings. For this reason, we expect that all radiative corrections to the tt+tt^ 7r°7r° 
scattering amplitude - due to transverse photons - vanish at threshold at order e^. We 
therefore neglect the transverse photons in the non-relativistic theory, while matching the 
relativistic and non-relativistic amplitudes at threshold. 
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APPENDIX C: CONTRIBUTION OF TRANSVERSE PHOTONS TO THE 

DECAY WIDTH 



We consider the role of transverse photons in the calculation of the decay width of the 
TT+TT^ atom. We evaluate their contribution for several typical diagrams and show that these 
do not contribute at order 5^/^. The procedure goes in several steps. 



1. As was mentioned in Section [111 the master equation ( |3.2(J| ) for the position of the 



bound-state pole is valid also in the case of a general non-relativistic Lagrangian. Expanding 
in a Taylor series around z = Eq gives 

_ (1'o|t(go)|*I-o) 

i-3l;(W£^o)l*o)"' • 

One may evaluate the denominator in this expression by retaining only leading contributions 
to t{Eo), given by strong bubbles with Coulomb ladders, 

(q|r(^)|p) = -ci - ^ p'/\z) cl + cl {gc{z)) + • • • , (C2) 

where p{z) = M^^o^z — 2Mt,o), and where the quantity {gc{z)) is defined in analogy with 
Eq. ( [4.6|) for the case of generic z. The explicit expression for this quantity is given by 

{gc(-)> = (gO-c(2)> + (g^c(;)> + {Sn-C(:)) , 

(go-c(j)) = ^ (A/,*(2M,+ - j))'/' , 

(8»-cW) = ^(*(2-'))-*(l)-V^). (C3) 

where ^(x) denotes the logarithmic derivative of Gamma-function. The quantity rj was 
defined after Eq. (|3.12|) , and A(/i) is given in Eq. ( p.28|) . Using the above expressions, it 
is easily seen that the width is modified at 0(5^^/^) in the presence of the denominator in 



Eq. (|CT|). Consequently, at the accuracy we are working, one may use 

r = -2Im2 = -2Im(*o|T-(^^o)|^o)- (C4) 

2. In general, the couplings in the non-relativistic Lagrangian are not real. Decay 
processes with an energy release at the hard scale contribute to the imaginary part of these 
constants. The only possible intermediate states in such diagrams are wy and tt^ + ri'j. 
Since the anomaly-induced decay into 7r° + 7 cannot proceed from the ground state due to 
C-invariance, the states with a minimal number of photons are 27 and vr" + 27. However, 
the decay width into two photons is of order 6^ ||l|,41| , and the decay width into vr" + 27 



starts, at least, at the same order in 5. Therefore, at order 5^/^ one may assume that all 
couplings in the non-relativistic Lagrangian are real, and that the Hamiltonian constructed 
from this Lagrangian is hermitian. 
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For a hermitian Hamiltonian, the operator r obeys the unitarity condition 

t{Eo) - t\Eo) = -27iit{E,) 6{Eo - Hq - He) t^Eo) , (C5) 

where the symbol 6 is defined as follows: in order to evaluate the right-hand side of Eq. (|C5|), 
one inserts a complete set of eigenstates (Hq + H^)]/?) = Ej^ljS), omitting the ground state 
of the of the bound tt+tt^ system. It is easy to see that the only allowed states are those 
containing either tt+tt^ + A^7, > or 27r° + A^7, = 2A; > scattering states (the decay 
into 2Tr^+ [odd number of photons] from the ground state is forbidden by C-invariance). 
The contribution from tt+tt" + A^7 vanishes due to lack of phase space. So we have 

r = T2nO + r2^o_|.27 + ^27t0+aj H (^al coupliugs) . (C6) 

It can be easily seen that the decay width into 2tc^ + 2'j starts at 0(5^^/^), and the decays 
into states containing four and more photons are even more suppressed. Consequently, 
T = T2^o + 0{5'). 

3. From the unitarity condition the following expression for the decay width into 2tt^ 
final state is readily obtained, 

r2.o = M^\T,^,^ko)\\ T2Ah) = I (^^o(p) +(p|T(Eo)|ko)o, (C7) 

where ko = |ko| = {M^o{Eo - 2iV4o)}i/2 ~ 0{6^/^) is the magnitude of the relative 3- 
momentum of the neutral pion pair in the final state, and the subscripts -|-, distinguish 
between charged tt+tt^ and neutral tt'^tt^ state vectors, respectively. 

The leading strong contribution in T2^o(A;o) starts at 0{5^^^). It is straightforwardly seen 
from Eqs. (|C^, ( |C3| ) and ( P?] ) that in order to evaluate width at 0(5^/^), it suffices to 
calculate T2„o(A;o) at 0(5^/^). Here we are interested in the contributions to this quantity 
due transverse photons. The diagrams that may potentially contribute in the lowest order 
in 6 are depicted in Fig. |^ - these are the self-energy (Figs. 0a,b) and vertex (Figs. ^,e) 
corrections to the lowest-order strong four-pion non-derivative vertex, with any number of 
Coulomb photons. In addition, there are the diagrams Fig. |^ that stem from the coun- 
terterm Lagrangian, Eq. (p^ ) - they are needed to renormalize the self-energy diagrams. 
These diagrams are the counterparts of diagrams depicted in Figs. ^,b, apart from the fact 
that in the latter there are no Coulomb ladders. One may consider the counterparts for 
other diagrams depicted in Fig. ^ and Fig. |^ which, however, are expected to be at least as 
suppressed in powers of 6 in the bound-state calculations, as the ones displayed in Fig. |^. 
Further, aiming to establish the power in S where these diagrams start to contribute, we 
may omit the Coulomb ladders: diagrams with Coulomb photons cannot be amplified with 
respect to the diagrams with no Coulomb photons. In order to prove this, we note that 
any additional Coulomb exchange in a given diagram adds an integration measure d^p, the 
Coulomb potential a\p — q|~^, and the energy denominator {z — ■ ■ ■ — 2M^+ — p"^ /M^'^. 
As the momenta scale like a or a^, depending on the topology of the diagram, it is seen 
straightforwardly that a diagram with an additional Coulomb photon, at worst, contributes 
to the same order in a as the original one. 

The contribution to the matrix element of r, coming from the diagram in Fig. |^ (with 
no Coulomb photons) is 
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[p|r(i^o)|ko)o = -2C2^^ = (^^^^j (4.)^/^r(l + i) + ■ 



^]o = 2M^+ + -^-Eo, (C8) 

where, in order to be consistent with the matching, we have used again the threshold ex- 
pansion. The quantity S(f2o,p^) in this expression is the self-energy part introduced in 
Eq. (p3|), calculated at the off-shell value Qq of the parameter Q - for this reason, one does 
not encounter an infrared divergence performing the limit d ^ 3. 

The contribution from the diagram Fig. ^a should be complemented with the counterterm 
contribution. In the Hamiltonian formalism, one has first to use the EOM in order to 
eliminate the time derivative in the counterterm Lagrangian (p6|). We find 



J f^f 



AH = -^-^ ld'^\j:{^iB ni){nono) + h.c\ + ■ ■ ■ . (C9) 



± 



The contribution from this Hamiltonian to the above matrix reads 

+ (p|AH|ko)o = I^P^- (CIO) 

As was expected, the contribution from the counterterm. Fig. 0c, cancels the UV divergence 
from the self-energy diagrams Fig. ^ and Fig. |^. 

Finally, folding this result with the ground-state wave function, and rescaling the integra- 
tion momentum as p 7p, it is seen that the total (self-energy + counterterm) contribution 
to T2,ro(^o) starts at 0(5^/^), and therefore can be neglected. Note that one should not be 
worried about the (spurious) UV divergence in the integral over the momentum p - these 



divergences cancel once all contributions in the given order in a are summed up |^| 
The contribution coming from the diagram in Fig. |^d is given by 



[p|^(^o)|ko 



e C2 f d^q p^q^ — (pq) 1 



|Fig.7d ^ 

° Ml^ j {2tiY |p-q|3 - 2M,+ - qVM,4 

1 



Eo - 2M,+ - |p - q| - pV(2M,+ ) - qV(2Af,4 
e^C2 r d^q^ p^q^ — (pq)^ 1 



M2+ J {27cY |p - qr ^0 - 2M^+ - q^M^, 



+ •••. (Cll) 



Again, the threshold expansion provides one with a homogeneous expression. After 
rescaling the momenta, it is immediately seen that the contribution to the quantity T2Tro{ko) 
starts at 0(5''/^) and therefore can be neglected. The contribution from the diagram ^ is 
identical to the one from Additional diagrams may be treated in an analogous manner 
- using power-counting arguments, we expect them to be even more suppressed. 

In summary, we have seen that - for a large class of diagrams - transverse photons 
contribute neither to the decay width at 0(5^/^), nor to the matching condition at threshold. 
Although we have not provided a mathematical proof, we believe that this result is true in 
general. For this reason, we completely neglect transverse photons in our non-relativistic 
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theory. The rest is then straightforward: one ehminates the Coulomb photons by using 
EOM, and retains only those terms in the Lagrangian that contribute to the decay width at 
0(5^/^). In this manner, one arrives at the Lagrangian displayed in Eq. ( |2.2| ). 



APPENDIX D: MATCHING PROCEDURES 

In this paper, we have matched the amplitudes in the relativistic and non-relativistic 
theories at physical space dimension d = 3. In this case, both amplitudes contain singular 
pieces that behave like |p|~^ and ln|p| in the vicinity of the threshold. In addition, there 
is the infrared-divergent Coulomb phase. The matching is performed at threshold, for the 
finite parts of the amplitudes that are obtained after removing the Coulomb phase, and 
subtracting the singular pieces. 

In the literature [0], there are examples of a different matching condition, where the 



matching is performed for the full amplitude at threshold at d ^ 3. The threshold singu- 
larities then, in general, transform into poles at d = 3 that cancel at a final stage. In this 
Appendix, we compare these matching conditions in two specific examples: we consider the 
two diagrams depicted in Fig. ^ and their non-relativistic counterparts. As we will see, the 
two matching conditions lead to exactly the same result. 

Let us start with the vertex correction depicted in Fig. ^a, that gives rise to the Coulomb 
phase, and to the |p|~^ singular behavior in the real part of the amplitude. In the non- 
relativistic theory, the corresponding vertex integral is given by Eqs. ( |2.24| ), ( |2.25| ) and 



( p.26| ). If one instead reverses the order of limiting procedures and calculates the same 



integral at p = 0, c? 7^ 3, one finds that the integral vanishes, K;(p, 2w{p)) = (the symbol 
"tilde" is used to distinguish the quantities calculated by using this particular sequence of 
limiting procedures). Both the Coulomb phase and the |p|~^ singularity disappear when 
this prescription is used. 

Let us now turn to the same diagram in the relativistic theory. The infrared-singular 
contribution at threshold is contained in the function G'_|__^(s) defined by 



rf^g 1 



(27r)«« q^{q^-2qp,){q^ + 2qp2) 



2 r 1-a . 

<^ In h '171 

sa I 1 + a 



1 A 



, ,27riln ^-4Li + 27r^ - In^ — — , (Dl) 

where pj = pi = M^+, s = {pi + p2f > AMl+, a = (1 - AMl+/sf'^, and 

Li(x) = Li2(l-x) = Y^y- 
This function was considered before in Ref. where the authors had introduced a photon 



mass to tame the infrared divergences. Here, we use dimensional regularization. [The 
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definition of Li2(x) in Ref. [0 has to be changed according to Eq. ( P2| )P|.] 



The infrared-divergent quantity \jr is defined by 

The rest of the diagram is infrared-finite at threshold. In the vicinity of the threshold, 

1 r 7raM^+ „ , a /, Ml 



If we reverse the sequence of limiting procedures, we find 

e'G^.,{s) ^ [l6naXjn + £ (^^ ^ + s) } • (D5) 

This procedure therefore again amounts to dropping the Coulomb phase and the singular 
piece that behaves like |p|~^. The matching condition is not altered since for the particular 
combination of Vc and G^^^ that appears in this condition, 

..1 A M! 



K - 4M^^G'+_^ = V; - 4M^^G+_-, = -AXiR - — ^\n-^ + 3j. (D6) 

Let us now consider the diagram depicted in Fig. |]b that leads to a logarithmic singularity 
at threshold. The corresponding loop integral in the non-relativistic theory is given by 
Eqs. ( |2.27| ) and ( |2.28| ). If one reverses the order of limiting procedures, one finds Bc{2w{p)) = 



0. 

In the relativistic theory, it suffices to consider the scalar integral 

Rc{P) = -e' J j0j,D^+{l)D^+{P - l)D^,{q)D^+{P - q)D{l - q) , (D7) 

with D^+{q) = (M^+ —q"^)^^, D{q) = (q^)^^- One need not consider diagrams with derivative 
couplings: those can be expressed through the integral Rc{P) and through integrals that are 
suppressed at threshold as compared to Rc{P), or are infrared-finite. 



The explicit expression for Rc{P) is given in Ref. |H6 



where D = 4 — 2e, and 

G,= + H,= 3F2(l,l + £,l + 2£;^ + £,2-£;^^), (D9) 

where, in difference with Ref. [^, P^ is defined in Minkowski space. 



^ We thank M. Knecht for correspondence on this point. 
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Rc{P) 



At D = 4, the integral near threshold in the CM frame (P = 0) behaves as 

2|p 



a 



1287rM2, 



2 In 



, o 21C(3) . 
+ ln2 + ^^-m 



0(|p|), P^ = {2^Ml^+^\Q). 

(DIO) 



If we first set p = and then consider the limit D ^ 4, we find 



Rc{P) 



a 



27r2 y 



, P^ = (2M^+,0), (Dll) 



with A(yu) given in ( |2.28| ). For the combination that appears in the matching condition, we 
have 



aMl 



57r 



27r2 ) 



■ (D12) 



Consequently, the matching condition remains unaffected by the interchange of the limiting 
procedures. 

To conclude, for these two diagrams we have checked that the result of matching is the 
same for the two prescriptions. We expect that this conclusion holds in the general case as 
well. 



APPENDIX E: THE MAPPING SU{?,) x SU{2,) SU{2) x SU{2) 

In order to perform the mapping SU{2) x SU{2) SU{3) x SU{3) for the constants 
ki and Ki, we evaluate the neutral pion mass and the amplitude tt+tt^ vr^vr" in the 
SU{3) X SU{3) framework, expand the result in powers of m/rris and compare it with its 
SU{2) X SU{2) analogue. 

From the expressions for the neutral pion mass, we find the relation 

lOkl + 10k; - 18kl + 9kl - lOkl - 10k; - 2k'j 

= 12KI + i2k; - i8k; + 9KI + iok; + iok; - i2k; - i2k; - iok; - iokIq . (ei) 

Matching the coefficients of s and rh in the tt+tt^ tt^tt^ amplitudes, we find 

lOkl - 8k; + 18kl - 9kl = 12KI - QK^ + 18K; - 9KI + IOK; - 8KI , (E2) 

and 

20kl - 16k; + 18k; - 9kl - 10k; + 26k; - 2kl^ + 36k; 
= 2AKI - 12KI + 18^3^ - 9KI + 207^5^ - 16KI - 12K^ + 2AKI - lOK^ 

+26K[o + 36ir[i - 144ZoL^ - 72ZqLI . (E3) 

Here, Zq is the SU (3) x SU (3) analogue of the coupling Z. In the order of the quark mass 
expansion considered here, we may identify Zq with Z. Combining the relations ([El|), ( [E2|) 
and (|E|), we obtain Eqs. and {W^ - 
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The coupling can be expressed p8| as a convolution of a QCD correlation function 
with the photon propagator, plus a contribution from the QED counterterms. We have 
checked that P{ki) in ( |4.26| ) is independent of the QCD scale /xq that must be introduced in 



the QCD Lagrangian after taking into account electromagnetic effects [|38| , |39 
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FIGURE CAPTIONS 



FIG. |l|. Examples of diagrams generated by the Lagrangian ( p.2|) at e = 0. Solid (dashed) 
lines correspond to charged (neutral) pions, crosses denote mass insertions, and the filled 
circle stands for a higher-order derivative vertex. 

FIG. |2|. Building blocks for the scattering amplitude, including Coulomb interactions at 
order e^. Dotted lines denote the exchange of a Coulomb photon. 

FIG. 1^. Diagrams contributing to the decay of the tt+tt^ atom into two photons in the 
relativistic theory. 

FIG. ^. Self-energy of the charged pion at O(e^). The twisted line denotes a transverse 
photon. The counterterm diagram b) stems from the Lagrangian (p6| ). 

FIG. ^. Radiative corrections to the scattering amplitude vr+vr" — > ir^n^, minimal couplings. 
The twisted lines denote transverse photons. Ellipses stand for any number of charged and 
neutral pion loops. Mass insertions are not displayed. 

FIG. 1^. The same as in Fig. |^, but with at least one vertex describing the coupling of 
a transverse photon to four pions (filled circles). Ellipses stand for any number of strong 
loops. 

FIG. 0. Self-energy and vertex corrections to the decay width of the tt+tt^ atom: (a,b) 
self-energy corrections, (c) contribution of the counterterm, (d,e) vertex corrections. The 
twisted lines stand for transverse photons. 

FIG. |8|. Diagrams that generate a singular behavior of the tt+tt^ tt^tt^ scattering ampli- 
tude at threshold in the relativistic theory: (a) vertex correction, (b) internal exchange of 
the photon. 



38 



FIGURES 






a 



FIG. 1. 





a 



FIG. 2. 



xAAA/> 



lAAA/^ 







FIG. 3. 



39 





FIG. 5. 

40 




41 



